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^ . Abstract 

Q ■ The fact that the Korteweg-de- Vries equation offers a good approximation of long-wave 

| solutions of small amplitude to the one-dimensional Gross-Pitaevskii equation was derived 

several years ago in the physical literature (see e.g. [17]). In this paper, we provide a rigorous 
proof of this fact, and compute a precise estimate for the error term. Our proof relies on the 
integrability of both the equations. In particular, we give a relation between the invariants 
of the two equations, which, we hope, is of independent interest. 



1 Introduction 



, In this paper, we consider the one-dimensional Gross-Pitaevskii equation 

►> . 

fC ■ idtV + dPV = ^(|^| 2 - 1) on K x R, (GP) 

t^J- ! which is a version of the defocusing cubic nonlinear Schrodinger equation and appears as a 

relevant model in various areas of physics: Bose-Einstein condensation, fluid mechanics (see e.g. 
[12, 18, 15, 6]), nonlinear optics (see e.g. [16]). 

qq ' We supplement this equation with the boundary condition at infinity 

|*(x,i)| -> 1, as Ixl -> +oo. (1) 

> ■ 

This boundary condition is suggested by the formal conservation of the energy (see (2) below), 
and by the use of the Gross-Pitaevskii equation as a physical model, e.g. for the modelling of 
"dark solitons" in nonlinear optics (see [16]). Note that boundary condition (1) ensures that 
(GP) has a truly nonlinear dynamics, contrary to the case of null condition at infinity where the 
dynamics is governed by dispersion and scattering. In particular, (GP) has nontrivial localized 
coherent structures called "solitons". 

At least on a formal level, the Gross-Pitaevskii equation is hamiltonian. The conserved 
Hamiltonian is a Ginzburg-Landau energy, namely 

= \ [ |d x *| 2 + - [ (1 - |*| 2 ) 2 = / e(W). (2) 
z Jr 4 Jr Jr 
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In this paper, we will only consider finite energy solutions to (GP). Similarly, as far as it might 
be defined, the momentum 

P(*) = I f (i&tf,*) (3) 
is formally conserved. Another quantity which is formally conserved by the flow is the mass 

m(*)^|(|f| 2 -l). (4) 

Equation (GP) is integrable by means of the inverse scattering method, and it has been 
formally analyzed within this framework in [19], and rigorously in [11]. The formalism of inverse 
scattering provides an infinite number of invariant functionals for the Gross-Pitaevskii equation 
and our proofs rely crucially on several of them. Concerning the Cauchy problem, it can be 
shown (see [20, 10]) that (GP) is locally well-posed in the spaces 



X k (R) = {u£ Ll oc (R, C), s.t. E(u) < +oo, and d x u £ H^ 1 

for any k > 1, and globally well-posed for k = 1. In the one-dimensional case considered here, it 
is also globally well-posed for k > 2. 

Theorem 1. Let k G N* and Vl/o £ A fc (R). Then, there exists a unique solution in 
C°(R,X fe (R)) to (GP) with initial data If *o belongs to X k+2 (R), then the map t ^ 

*(•,£) belongs to C 1 (R,X k (R)) and C°(R, X k+2 (R)). Moreover, the flow map ^ ^ ®{;T) is 
continuous on X k (R) for any fixed T € R. 

Furthermore, the energy is conserved along the flow, as well as the momentum, at least under 
suitable assumptions (see e.g. [3]). On the other hand, the rigorous derivation of conservation 
of mass raises some difficulties. 

If <I> does not vanish, one may write 

* = y^expit/?. 

This leads to the hydrodynamic form of the equation, with v = 2d x f, 

( d t p + d x (pv) = 0, 

\d, P \ 2 \ (5) 



p(d t v + v.d x v) + d x (p 2 ) = pd x {^ - J^P) . 



If one neglects the right-hand side of the second equation, which is often referred to as the 
quantum pressure, system (5) yields the Euler equation for a compressible fluid, with pressure law 
p(p) = p 2 . Since the right-hand side of (5) contains third order derivatives, this approximation 
is only relevant in the long- wave limit. A rigorous derivation of this asymptotics was derived by 
Grenier in [14] for different conditions at infinity. 

Recall that linearizing the compressible Euler equation with pressure p(p) = p 2 , around the 
constant solution p = 1 and v = 0, one obtains the system 

<9 t p + <9 x O = 0, () 
dtX> + 2d x p = 0, w 

which is equivalent to the wave equation with speed c s given by 

c 2 s = 2. 
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This speed is referred as the sound speed for the Gross-Pitaevskii equation. In this setting, the 
wave equation (6) appears as an approximation of the Gross-Pitaevskii equation. As mentioned 
above, this amounts however to neglect the quantum pressure, coming from the dispersive prop- 
erties of the Schrodinger equation, as well as to restrict ourselves to small long-wave data, so 
that the wave equation approximates the Euler equation. Rigorous mathematical evidence of 
this fact is provided in [1]. 

In order to specify the nature of the perturbation as well as of the long-wave asymptotics, 
we introduce a small parameter < e < 1 and set 

f p(x,t) = 1 + -^a £ (ex,et), 
| v(x,i) = ev £ (ex,et), 

so that system (5) translates into 

d t a £ + \f2d x v £ = -ed x (a £ v £ ), 
' d t v e + V2d x a E = e(- VE . dx v 6 + 2d x ( ) • (?) 

Specifying a result of [1] in dimension one, we are led to 

Theorem 2 ([1]). Let s > 2. There exists some positive constant K(s) such that, given any 
initial datum (a £ ,v £ ) G fF +1 (R) x H S (R) verifying 

K(s)e\\(a £ ,v £ )\\ Hs+ i (R)xH s (R) < 1, 

there exists some real number 

1 

T £ > 



K(s)e 2 \\(a £ ,v^)\\ H s+l mxH s m , 
such that system (7) has a unique solution (a £ ,v £ ) £ C°([0, eT £ ], H S+1 (R) x H S (W)) satisfying 

\\(a £ (-,et),v £ (-,et))\\ H s+i {R)><H s {R) < K(s)\\(a° £ ,v° £ )\\ H s+i m><H s m , and i < p(-,t) < 2, 

for any t £ [0,T £ ]. Let (a, o) denote the solution of the free-wave equation 



(8) 



d t a + V2d x v = 0, 
dtt> + V2d x a = 0, 

with initial datum (a £ ,v £ ), then for any < t < T £ , we have 

\\(a £ ,v £ )(-,et) - (a,D)(-,et)|| jf/s -2 (R)x ^ s -2 (R) 

< K(s)(e 2 t\\(a £ ,v°)\\ 2 Hs+1(R)xHs(R) + e 3 t||(a°, ^)||^ s +i (R)xHs(R) ) . 

Remark 1. Notice that the bounds on K(s) provided by the proof of Theorem 2 in [1] blow up 
as s tends to +oo. An interesting question is therefore to determine whether the constant K{s) 
may be bounded independently of s. In particular, it would be of interest to extend the result 
to the limiting case s = +oo. 

The purpose of the present paper is to consider even smaller perturbations of the constant 
one, and to characterize the deviation from the wave equation on larger time scales. Our initial 
data has the form 

r P (x,o) = i-^iv £ o ( e x), 

\ u(x ,0) = ^°( £ x), 
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where N® and W® are uniformly bounded in some Sobolev space H S (R) for sufficiently large s. 
Applying Theorem 2 to such data, that is for a £ = — ^pAf? and v £ = W £ , yields uniform 
bounds on a time scale T £ = 0(e~ 3 ). More precisely, setting 

6 6\/2 

n e (ex, et) = -=a e (sx, et), and w e (ex, et) = v e (ex.,et), 

eV2 £ 

it follows for such initial data from Theorem 2 that we have 

Proposition 1. Assume s>2 and Ke 2 \\(N®, W £ )\\ HS+ i(u) xH3 (^ < 1. Let (n, tt>) denote the 
solution of the free wave equation 



(9) 



(10) 



dt(y/2n) -d x tt) = 0, 

a t «-2a x (V2n) =0, 

initial datum (AT?, W° ). Tnen, /or any < t < T £ , we have 

\\{n £ ,w £ )(-,et) - (n,tt>)(-,et)|| H .-2 (R ) xH .-2 (R ) 
< Keh(\\(N°,W?)\\ HS+ i (R)xH s {R) + W"°)|&. +1(R)xffs(I 

where 

T = I 

£ Ke*\\{N2,W*)\\ H s +imxH s m 

In particular, if N £ and W £ are required to be uniformly bounded in .fP +1 (IR) x H S (M), then 
in view of (10), the wave equation provides a good approximation on time scales of order o(e~ 3 ). 
This approximation ceases to be valid for times of order C(e~ 3 ) as the subsequent analysis will 
show. 

The general solution to (9) may be written as 

(n,tt>) = (n + ,tt> + ) + (rT,rtr), 

where the functions (n ± ,rD ± ) are solutions to (9) given by the d'Alembert formulae, 

(n+(x,t),ro + (x,t)) = (N + (x-V2t),W + (x-V2t)), 
(n _ (x,t),tt)-(x,t)) = (N-(x + V2t),W-(x + V2t)), 

where the profiles and are real-valued functions on R. Solutions may therefore be split 
into right and left going waves of speed \f2~. Since the functions (n^ 1 , ro ± ) are solutions to (9), it 
follows that 

(27V+ + W + ) x = 0, and (2N~ - W~) x = 0, (11) 
so that, if the functions decay to zero at infinity, then 

2JV ± = tW ± = — — £ -. (12) 

At this stage, it is worthwhile to notice that the Gross-Pitaevskii equation, as well as the wave 
equation, is invariant under the symmetry x — > —x. 

It remains to derive the appropriate approximation for time scales of order 0(e~ 3 ). On a 
formal level, this was performed in [17]. We wish to give here a rigorous proof of that approxi- 
mation. In view of the previous discussion, and following the approach of [17], we introduce the 
slow variables 

e 3 

x = e(x + y/2t), and r = — -=t. (13) 

2\/2 
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The definition of the new variable x corresponds to a reference frame travelling to the left with 
speed y/2 in the original variables (x, t). In this frame, the wave (n _ ,tu _ ), originally travelling 
to the left, is now stationary, whereas the wave (n + , ro + ) travelling to the right now has a speed 
equal to 8e~ 2 . This change of variable is therefore particularly appropriate for the study of 
waves travelling to the left. This will lead us to impose some additional assumptions which will 
imply the smallness of N + and W + . Notice that the change of frame breaks the symmetry of 
the original equations. 

In view of (13), we then define the rescaled functions N £ and B e as follows 



6 , , 6 fx 4r 2>/2t 
N £ (x,t) = ^??(x,t) = - — 



r2''V-'"> £ 2''V £ £ 3' £ 3 

„ , , 6^2 , . 6^2 fx At 2y/2r- 

where \& = gexpiip and r/ = l — g 2 = l — \^\ 2 - 
Our main theorem is 



(14) 



Theorem 3. Let e > be given and assume that the initial data ^o(') = ^( - >0) belongs to 
V 4 (R) and satisfies the assumption 

\\N°\\ H s {R) + e||^iV £ || i2(R) + \\d x e° £ \\ H s m < K . (15) 

Let N £ and M. £ denote the solutions to the Korteweg-de Vries equation 

d T N + dlN + Nd x N = (KdV) 

with initial data N®, respectively d x O®. There exists positive constants eq and K±, depending 
possibly only on Kq such that, if e < £q, we have for any r £ R, 

We(;r) - N £ (;t)\\ l2{r) + \\M £ (-,r) - d x G £ (; r)|| L2(K) 
<K 1 (e+ \\N° - d x e £ \\ H s m ) exp^xlrl). 

Theorem 3 yields a convergence result to the (KdV) equation for appropriate initial data. 
Since the norms involved in (19) are translation invariant, the (KdV) approximation can only 
be relevant if the waves travelling to the right are negligible. In view of our previous discussion, 
this is precisely the role of the term \\N® — d x G £ \\fj3m) in the right-hand side of (19). Indeed, 
in the setting of Theorem 3, the right going waves iV + and W + are given by 

2N + = -W + = iV £ ° - d x @° £ . 

If the term \\N® — d x Q £ \\ H s^ is small, then the (KdV) approximation is valid on a time interval 
(in the original time variable) t 6 [0, S £ ] with 



S F = cm min 



log(e)| | log(||iV0 - ^eg||^ 3( 



E 6 E 6 



In particular, if \\N® — d x @ £ \\ H 3^ < Ce a , with a > 0, then the approximation is valid on a 
time interval t € [0,S £ ] with S £ = o(e~ 3 | log(e)|). Moreover, if \\N® — d x Q^\\ H 3^ is of order 
0(e), then the approximation error remains of order 0(e) on a time interval t G [0, S e ] with 
S £ = 0(e-1). 



5 



Remark 2. We also show in the course of our proofs (see Proposition 2 below) that, under the 
assumptions of Theorem 3, the i? 3 -norms of N £ and d x Q £ remain uniformly bounded in time. 
Since the same property holds for the solutions M £ and Ai £ , it follows by interpolation that the 
difference of the two solutions may also be computed in terms of .fP-norm as 

\\K(;r) - N £ (;r)\\ H s {m + ||M e (-,r) - d x @ £ (; r)\\ H s {R) 
< K(e+\\N°-d x e° £ \\ HHR) ) a{s) e X p(a( S )K 1 |r|). 

for any < s < 3 and any r G R, where a(s) = 1 — |, and where the constant K depends 
possibly on Kq and s. 

Remark 3. As a matter of fact, we believe that for any s > 0, the following inequality holds 

\\K(;r) -N £ (;t)\\ H s {r) < K(s)(e 2 + ||iV e - d x @° £ \\ Hs+ s {R) ) ex P (^i|r|), (20) 

for any r G R. To prove inequality (20) along the lines of the proof of Theorem 3 would require a 
more general treatment of the invariants of the Gross-Pitaevskii equation, whereas in this paper, 
we have only handled the lower order ones (at the cost of sometimes tedious computations). In 
a forthcoming paper [4], we make use of a different strategy avoiding invariants but at the cost 
of a higher loss of derivatives. Here also, as in Remark 1, it would be of interest to prove a result 
in H°°(R). 

The functions N £ and d x Q £ are rigidly constrained one to the other as shown by the following 

Theorem 4. Let ^ be a solution to (GP) in C°(R,H 4 (R)) with initial data Assume that 
(15) holds. Then, there exists some positive constant K , which does not depend on e nor t, such 
that 

||iV £ (-,r)±^e £ (-,r)|| L2(M) < ||7V0±^e°|| i2(M) +^ £ 2 (l + |r|), (21) 

for any r G R. 

The approximation errors provided by Theorem 3 and 4 diverge as time increases. Concerning 
the weaker notion of consistency, we have the following result whose peculiarity is that the bounds 
are independent of time. 

Theorem 5. Let ^ be a solution to (GP) in C°(R, if 4 (R)) with initial data Assume that 
(15) holds. Then, there exists some positive constant K , which does not depend on e nor r, such 
that 

\\d T u £ + dlu £ + u £ d x u £ \\ L2m < K(e + ||iv e - d x e° £ \\ Hm) ), (22) 

for any r <= R, where U £ = N ^+^ e ^ . 

The relevance of the function U £ will be discussed below. 

A typical example where the assumptions of Theorem 3 apply is provided by travelling wave 
solutions to (GP), i.e. solutions of the form ^(x, t) = v c (~x + ct), where the profile v c is a 
complex-valued function defined on R satisfying a simple ordinary differential equation which 
may be integrated explicitly. Solutions then do exist for any value of the speed c in the interval 
[0,\/2). Next, we choose the wave-length parameter to be e = \Jl — c 2 , and take as initial data 
^ e the corresponding wave v c . We consider the rescaled function 
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where rj c = 1 — \v c \ 2 . The explicit integration of the travelling wave equation for v c leads to the 
formula 

y £ (x) = y(x) = — „ . . . 
£K ' K ' ch 2 (f) 

The function v is the classical soliton to the Korteweg-de Vries equation, which is moved by the 
(KdV) flow with constant speed equal to 1, so that 

N £ (x, t) = v(x - r). 

On the other hand, we deduce from (14) that = v, so that 



4 / I e 2 ^ 



N £ (x,t) =y(x- -^(l - yi 

Therefore, we have for any rGl, 

||A4(-,r)-iV £ (.,r)|| L2(M) = 0(e 2 r). 
Concerning the phase (p c of i> c , we consider the scale change 

o 6^2 fx\ 

so that, in view of [13], 

eV ' V 2 l-fv(a 

and hence, 



\\N»-d x ®l\\ HHm = 0{e 2 ). 

This may suggest that the e error in inequality (19) is not optimal. As a matter of fact, we 
believe that the optimal error term would be of order e 2 (as mentioned in formula (20)). A proof 
of this claim would require to have higher order bounds on N e and d x 6 e . 

We next present some ideas in the proofs. We infer from (GP) the equations for N £ and £ , 
namely 

d x N £ - d 2 e £ + £ - {^d T N £ + i N £ d 2 x e £ + hd x N £ d x e £ ) = o, (23) 

and 

9e -N +-(-d@ I &1 I l ( 5 Q)A I g4 (MeI -q (24) 

6 £ ^ 6 6 ' 

The leading order in this expansion is provided by N £ — d x Q £ and its spatial derivative, so that 
an important step is to keep control on this term. In view of (23) and (24) and d'Alembert 
decomposition (12), we are led to introduce the new variables U £ and V £ defined by 

a^ £ + d x e £ n £ - d x & £ 

U £ = , and V £ = , 

and compute the relevant equations for U £ and V £ , 

d T U £ + dlU £ + U £ d x U £ = -d x V £ + U x (U £ V £ ) + l -d x (V 2 ) - e 2 R £ , (25) 

and 

d T V £ + \d x V e = d 3 x U £ + 8lV £ + \d x {V 2 ) - \d x {U £ ) 2 - \d x {U £ V £ ) + e 2 R £ , (26) 
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where the remainder term R £ is given by the formula 

= N £ dlN £ (d x N £ )(d 2 N £ ) e 2 (d x N £ f 
£ 6(l-fiV £ ) 3(l-fiV £ )2 36(l-f7V £ )3- 

The left-hand side of equation (25) corresponds to the (KdV) operator applied to U £ : a major 
step in the proof is therefore to establish that the right-hand side is small in suitable norms. This 
amounts in particular, as already mentioned, to show that V £ , which is assumed to be small at 
time r = remains small, and that U £ , which is assumed to be bounded at time r = 0, remains 
bounded in appropriate Sobolev norm. To establish these estimates, we rely among other things 
on several conservation laws which are provided by the integrability of the one-dimensional (GP) 
equation. To illustrate the argument, we next present it for the L 2 -norm, where we only need 
to invoke the conservation of energy and momentum. 

In the rescaled setting, the Ginzburg-Landau energy may be written as 

«4(/J(w !+ *')4/i^4« e '> , ))4* 9 - 1 ' (28) 

so that assumption (15) implies that 

S 1 (N^,e° e )<K . (29) 

On the other hand, when the energy E(*&) is sufficiently small, which is the case at the limit 
e — > 0, we may assume that 

^<|*|<2, 

which may be translated as 

\ < 1 - jN £ < 4, (30) 

so that the rescaled energy £ i satisfies 

I ((d x Q £ ) 2 + N 2 ) < KS^Ne, 8 £ ), (31) 

JR v 7 

where K is some universal constant. Similarly, the momentum may be written as 

p W = 11^ = ^1*^ s £w,a.>. (32) 

Next, we compute 

If r 2 f / 8 N 2 1 

£ 1 (n £ ,q £ )-V2v 1 (n £ ,q £ ) = - / {N e - d x e £ f + - / - * e - -N £ (d x e £ 

° Jr ° Jr\1- ^-N £ 3 

so that 

£i(N°,e° £ ) - v^i(JV?,e°)| < K . (33) 
Moreover, by the Sobolev embedding theorem and the inequality 2ab < a 2 + b 2 , 

£i(AT £ ,G £ ) - V2Pi(iV £ ,9 £ ) >K~ l U V £ 2 + e 2 [ (d x N £ f) - Ke 2 ( f (d x e £ ) 2 \ , (34) 



where K refers to some universal constant. By conservation, we then have 

^(£i(N £ ,e £ ))=0, and ± 



d (fi(iV e , 9 £ )) = 0, and ^{Pi(N £ , 9 £ )) = 0. (35) 
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Invoking (29) and (31), we are led to 

ll^(-^)ll! 2 (M) + ll^e £ (-,r)||| 2(K) <ir , 

for any r£l. In turn, using (33), (34) and (35) yields 

\\V £ (-,r)\\l 2m <K{\\V £ (0)\\l2 m +e 2 ). 

It turns out that the other conservation laws for the Gross-Pitaevskii equation involve quan- 
tities which behave as higher order energies and others which behave as higher order momenta. 
We denote £k(N £ , Q £ ) and T , k(N £ , Q £ ), respectively these quantities (precise expressions are pro- 
vided in Section 3). Using these invariants, we may perform a similar argument to control higher 
Sobolev norms. This gives 

Proposition 2. Let ^ be a solution to (GP) in C°(R, H 4 (R)) with initial data Assume 
that (15) holds. Then, there exists some positive constant K , which does not depend on e nor 
t, such that 

\\N £ (;r)\\ H s {R) +e\\d 4 x N £ (;r)\\ LHR) + \\d x e £ (;r)\\ H s {R) < K, (36) 

and 

\\N £ (.,T)±d x e £ (;r)\\ HHm <K{\\N°±d x G° £ \\HHR)+£), ( 37 ) 

for any TGl. 

The proof of Theorem 5 follows directly from Proposition 2. Using a standard energy method 
applied to the system (23) and (24) and taking advantage of the fact that the left-hand side of 
equation (24) is a transport operator with speed p-, we obtain Theorem 4. Finally, the proof of 
Theorem 3 follows again from an energy method applied to the difference W £ = N £ — N £ (and 
the equivalent for d x Q £ ). 

Remark 4. It is worthwhile to stress that in the course of proving Proposition 2, we have been 
led to prove a number of facts which, we believe, are of independent interest, and represent 
actually the bulk contribution of our paper. First, we have given expressions of the invariant 
quantities and proved that they are well-defined on the spaces X k (R): their expressions are not a 
straightforward consequence of the inductive formulae for the conservation laws provided by the 
inverse scattering method of [19]. Indeed, various renormalizations have to be applied to give a 
sound mathematical meaning to the expressions. Moreover, we have rigorously established that 
these quantities are conserved by the (GP) flow in the appropriate functional spaces. 

In a related direction, we have highlighted a strong and somewhat striking relationship be- 
tween the (GP) invariants and the (KdV) invariants. More precisely, we have shown that, for 
any 1 < k < 4 and for any functions in the appropriate spaces, 

£ k (N,d x Q) - V2V k (N,d x e) = g fc W ( Ar ~^ 9 ) + 0(e 2 ), 

where Ej? dV refers to the (KdV) invariants (for more precise statements, see Proposition 4.2). 
In particular, the (GP) invariants £ k and Vk, as well as the (KdV) invariants, provide control 
on the i? fc -norms. 

Remark 5. It would be of interest to investigate further the relationships between (GP) and 
(KdV), in particular at the level of the spectral problems associated to the corresponding inverse 
scattering methods. Indeed, recall that (KdV) can be resolved using scattering and inverse 
scattering methods for the linear Schrodinger equation 
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whereas (GP) is known to be tractable using the scattering and inverse scattering methods for 
the Dirac operator 



'1 + V3 \ fd x *i\ ,/0 /$i 

- 2 / "I* o J U 2 



Besides, it is known that the Schrodinger equation is a nonrelativistic limit of the Dirac equa- 
tion. Kutznetsov and Zakharov [17] suggest that this correspondence can be carried out in the 
asymptotic limit considered here. Notice however that rigorous scattering and inverse scattering 
methods require decay and regularity assumptions on the data (see e.g. [11] where the datum 
is required to decay at least clS I X I ; clS well as its first three derivatives). 

Let us emphasize again that our paper focuses on the left going waves. Our proof requires to 
impose conditions on the initial data to ensure that the right going wave is small. An interesting 
problem is to remove this assumption, i.e. to consider simultaneously both left and right going 
waves, and to study their interaction. We hope to handle this problem in a forthcoming paper, 
as well as the already mentioned optimal bounds. 

The paper is organized as follows. The next section is devoted to properties of the Cauchy 
problem. In Section 3, we compute the invariants of the (GP) flow needed for our proofs, 
and show that they are conserved. In Section 4, we recast these invariants in the asymptotics 
considered here, and show the convergence to the (KdV) invariants. In Section 5, we give the 
proofs to Proposition 2 and Theorem 5. Finally, in Section 6, we present the energy methods 
which yield the proofs to Theorems 3 and 4. 

While completing this work, we learned that D. Chiron and F. Rousset [5] were obtaining 
at the same time several results which are related to our analysis of the (KdV) limit, and also 
treated the higher dimensional case. 
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2 Global well-posedness for the Gross-Pitaevskii equation 

The purpose of this section is to present the proof of Theorem 1. It is presumably well-known 
to the experts, but we did not find it stated in the literature, and therefore we provide a proof 
here for the sake of completeness. 

Notice that Gallo [7] already established the local well-posedness of (GP) in the spaces A fc (M) 
for any k > 1 (see also [20, 10]). More precisely, we have 

Theorem 2.1 ([7, 10]). Let k > 2. Given any function S X k (R N ), consider the unique 
solution to (GP) in C°(R,X 1 (R)) with initial data * - Then, there exist (T_,T+) G 

(0,+cx)] 2 such that the map t ^ f(.,t) belongs to C°((-T_, T+), X k (R)). Moreover, either T + 
is equal to +oo, respectively T_ = +oo ; or 

\\d x V(;t)\\ Hk -i {R) +oo, as t T+ (resp. t -T_). (2.1) 
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7/ tf belongs to X k+2 (R), then the map t ^ *(■,£) belongs to C 1 ((-T^,T + ),X k (R)) and 
C°((-T_,T + ),X fc+2 (R)). Moreover, me /low map * ^ is continuous on X k (R) for 

any fixed —71 < T < T + . 

In view of Theorem 2.1, the proof of Theorem 1 reduces to establish that the i7 fc_1 -norm of 
the function d x *& cannot blow-up in finite time. In [7, 9], it is proved that the linear Schrodinger 
propagator S(t) maps X k (R) into X k (R), so that we may invoke the Duhamel formula 

*(•, t) = S(t)V - /' 5(/ - *)*(■, S) (1 - |*(-, S)\ 2 )ds, 

Jo 

to estimate the 7f fc_1 -norm of the function d x ^ by 



l^x*(-,t)||iffc-i(R) < ||5 x * ||iffc- 



+ 



f \\d^(-,s){l - \y(.,s)\ 2 ))\\ Hk - lm ds 
•J 



(2.2) 



To estimate the second term on the left-hand side, we invoke the following tame estimates. 

Lemma 2.1. Let k > 1 and (^1,^2) G X fc (R) 2 . Given emu 1 < j < k, there exists some 
constant K(j,k), depending only on j and k, such that 



< ^^(ii^iIIl-wII^II^w + II^||l~ ( «)RVi||l 2W )- 

We postpone the proof of Lemma 2.1 and first complete the proof of Theorem 1. 
Proof of Theorem 1. In view of (2.3), inequality (2.2) yields 



(2.3) 



\d^(-,t)\\ H k-i m < \\d x y \\ H k-i m +K(k) 



I (l+ll*(-,s)|| 2 oo (K) )||5 x *(-,s)|| fffc -i (K) ds 
J 



, (2-4) 



where K(k) is some constant depending only on k. Notice that *o belongs to X 1 (R), so that 
in view of the conservation of energy proved in [9] (see Theorem 3.1 below), we have 

E(*(;t))=E(* ). 

Next, given any function tp G X X (IR), there exists some universal positive constant K such that 



< 



(2.5) 



In particular, it follows from (2.5) that ||*(s)||loo( K ) < K (l + E(^ )) 2 , so that by (2.4), we are 
led to 



\d x *(-,t)\\ H k-i m <K(k,*o)[l + 



J \\d x ^{-,s)\\ H k-i {R) ds ^. 



where K (k, *o) is some constant only depending on k, E(^>o) and ||3 x VE r o||.fffc-i( 
we have by integration 

||<9 X *(-, t)\\ Hk - 1{R) < K(k, * ) exp (K(k, * )|t|), 
and it follows, going back to (2.1), that 

T_ = T+ = +00, 

which completes the proof. 



Therefore, 



□ 
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We now provide the proof of Lemma 2.1. 
Proof of Lemma 2.1. We introduce some cut-off function \ £ C°°(R, [0, 1]) such that 

X = 1 on (-1, 1), and X = on R \ (-2, 2), (2.6) 

and set 

X*(x) = x(-|), VxGR ' ( 2 - 7 ) 
for any R > 1. Using standard tame estimates, we have 

||^(XflV'iXfl^2)|| L 2( R ) 

< ^(i,^)(llXi?V ; i||L°°(iR)l|5x(Xi?V ; 2)||L2 (M) + ||x^2||l°o(r)||^(x^i)IIl2(r)) (2.8) 

< K(j, fc)(||^l|| L <»(R)||^(XJi^2)llL3(R) + ll^2|U~(R)ll^(XJi^l)lll^(R))- 

We now claim that 

R(Xfl^)IUacR) - II^H^R), as i? - +oo. (2.9) 

for any function ^ £ X fc (M) and any 1 < j < k. As a matter of fact, by the Leibniz formula, we 
have 

4M) = E Cf^x^r^- (2-10) 

m=l 

We next deduce from the dominated convergence theorem that 

XRd& ^ in £ 2 W, as -► +oo, 
whereas, when m > 1, we similarly have using (2.6) and (2.7), 

^^=Tll^" m ^lli«»(R) - 0, as i? ^ +oo. 
Hence, in view of (2.10), we are led to 

di(XR^) -» m L 2 (M), as R +oo, 

which ends the proof of claim (2.9). Combining (2.8) with (2.9), and noticing that (2.9) remains 
valid replacing xr by x%, we obtain (2.3) at the limit R — > +oo. This concludes the proof of 
Lemma 2.1. □ 

3 Invariants of the Gross-Pitaevskii equations 
3.1 Formal derivation of the invariants 

In [19], Shabat and Zakharov established that the one-dimensional Gross-Pitaevskii equation 
is integrable, and admits an infinite number of conservation laws f n (^>), leading to an infinite 
family of invariants / n (\I/). Set 

/!(*) = 4|*| 2 . (3.1) 
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and let 

/ n+1 (*) =^ x (^p) /,-(*)/„_,-(*). (3.2) 

Using the inverse scattering method, it is shown formally in [19] that the functions f n (^) ar e 
conservation laws for (GP), so that the related integral quantities I n (^) defined by 

4W= / (/„(*)(x)-/ n (*)(oo))dx, (3.3) 

are invariants for (GP). Here, the notation / n (\I/)(oo) stands for the limit at infinity of the map 
f n (^) assuming that 

Vl/(x) — > 1, as |x| — > +oo, and d^(x) — > 0, as |x| — > +oo, 

for any fc£N*. The first five conservation laws are computed in [19], namely (3.1) and 

/ 2 (tt) = -^VdxV, (3.4) 
/3W = -^* + i|*| 4 , (3.5) 
/ 4 («r) = --*I><9 X * + |^| 2 ^a x ^ + j\y\ 2 ^dj>, (3.6) 

/ 5 (tf) = -^<9** + ^|*| 2 *a 2 * + ^| 2 ^a 2 ^ + ^| 2 |a x ^| 2 + ^) 2 (d x ^) 2 - i|^| 6 . (3.7) 

The purpose of this section is to give a rigorous meaning to these quantities, to prove that 
they are conserved, and to extend the explicit list of invariants. As a matter of fact, these 
invariants enter directly in our analysis of the transonic limit. 

The first step is to compute the additional conservation laws using formula (3.2). Notice first 
that formula (3.2) is singular at the points where ^ vanishes. A first task is therefore to show 
that (3.2) can be used to define the functionals f n (&) even in the case the function ip vanishes 
somewhere. To remove the singularity in (3.2), we check by induction that the function f n (^) 
may be written as 

/„(*) = * x .F n (¥), (3.8) 
where the map T n is inductively defined by 

= (3.9) 

and 

71—1 

.F„+i(tf) = d x F n (V) + *2^-(tf).F n _ j (tf). (3.10) 

i=i 

In particular, the map J- n (^) is a polynomial functional of the functions • • • , <9 X _2 \E', 

Qn-2\p anc [ d™ -1 ^ which is defined without additional assumptions on \£. This leads to explicit 
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expressions of fe{^), .M 1 ^), /s(^) an d fg(^), which are given by 
/ 6 (tf) = - ^tfSgtf + 2|* \ 2 Hd^ + ^| 2 ^d 3 1> + 2\^\ 2 d x ^d 2 ^ + 3\$>\ 2 djrd^ 

+ ^W 2 ^x^^ + y |d x #| 2 *d x * - ^| 4 ^a x * - 2|^| 4 ^5 X *, 

_ in to 

+ 7(^) 2 9 X ^5 X * + 5|^| 2 |a 2 *| 2 + — (a x *) 2 *a 2 ^ + — (*) 2 (s 2 ^) 2 + i3|a x ^| 2 ^a 2 ^ 

U Q _ OK 

- |^|%5 2 ^ - y |^|%a 2 ^ - j^| 2 0f ) 2 (d x ^) 2 - 8|^| 4 |a x ^| 2 - y |*| 2 (i>) 2 (d x ^) 2 

/ 8 (*) = _ ^mlv + ^\ 2 ^d^ + 3|^| 2 ^a^ + 3|^| 2 a x ^a x ^+ y|^| 2 a x ^a x ^ 
+ io(^) 2 9 x ^9 x ^ + y|^| 2 9 2 ^9 x "^+ ^(9 x ^) 2 "^a x ^ + io|^| 2 a 2 ^a x ^ 

_ 55 71 — 5 

+ 17(^) 2 9 2 ^^ + 25\d x ^\ 2 ^d x ^ + —\d^\ 2 ^d x ^ + — (d 2 V) 2 Vd^ - | 4 ¥5gtt 

5 _ 5Q 75 

- 6|^| 4 *a x * - ^| 2 (^) 2 a x ^9 2 ^ - y |^| 4 d x ^s 2 ^ - y |^| 4 d x ^ d 2 ^ 
-27|*| 2 (^) 2 a x ^a 2 ^- ^|^| 2 |9 X ^| 2 ^9 X ^- ^\y\ 2 \d x y\ 2 vd x y - y($) 3 (9 s $) 3 

29 a — a — 

+ —\m\Hdy^ + 4|*| 6 *d x *, 

and 

/ 9 (*) = _ + J|^| 2 ^* + ^|*| 2 *a^ + ^| 2 a x ^a x ^+ y |^| 2 a x *a x ^ 

97 21 41 3^ 

+ y W 2 a x ^a x ^ + — 1^| 2 9 2 ^9 4 ^ + — (9 x ^) 2 *a 4 ^ + y |*| 2 d 2 *d 4 * 
+ ^ (*) 2 s 2 *a 4 * + y |c> x *| 2 ^<9 4 ^ + y |^| 2 |9 X *| 2 + y^a x ^a 2 ^9 x ^ 
+ y (*) 2 (s x *) 2 + ^a^s^ag* + Ip*^ 2 ^* + H|s 2 *| 2 ^ 2 * 

q _ q 15 70 

- ^| 4 f d 4 f - — 1*| 4 ^5 4 ^ - — |^| 2 (^) 2 9 x ^ 3 * - — |f | 4 <9 X ^<9 3 ^ 

- 36|17| 4 d x ^d x ^ - 49|^| 2 (*) 2 d x ^d 3 ^ - ^| 2 (^) 2 (9 2 *) 2 - ^-\^\ 4 \d 2 V\ 2 

1«5 1 QQ 

_ ^-\^\ 2 \d^\ 2 ^d^ - 66|^| 2 ^(d x ^) 2 d 2 ^ - ^|^| 2 (^) 2 (5 2 ^) 2 

- 29|^| 2 ^(a x ^) 2 a 2 ^ - ^|*| 2 |3 X *| 2 ^ 2 * - ^t^f(d x *fd 2 * - ^|*| 2 |d x *| 4 

- ^\d x ^\ 2 m 2 (d^) 2 + ^\^d 2 * + y \*\ e m 2 ^ + 7 -±\*\\M,)\d x W) 2 

+ ^i*i 6 i9 x *i 2 + wow 2 - ^i*r 

The second step is to provide explicit expressions of the invariants I n (^) associated to each 
conservation law f n {^) for an arbitrary function ^ in the appropriate space. This raises 

some serious difficulties since the integrands are not in general integrable when vl/ belongs to 
For instance, according to definition (3.3), the invariants h(^), h{^) and I^W) should 
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be given by 

im = \! {l-\*\ 2 )j2(*) = -\ / and/ 3 (*) = -J !*%* + \ I (i^i 4 -i). 

z jr z jr z Jr 4 jr 

(3.11) 

For an arbitrary function ^ in none of the above integrands belong to L 1 (M). Some 

quantities like ^fd 2 ^ can be handled using integration by parts. This is not possible for 1 — |\I>| 2 
or |^| 4 — 1, which do not involve derivatives. Even the quantity ^>d x ^ cannot be immediately 
treated by integration by parts. In particular, the renormalization process as used in formula 
(3.3) is not sufficient to give a sense to the invariants I n (^) in the spaces X k< 



When n = 2m + 1 is an odd number, a simple way to remove this difficulty is to introduce 
linear combinations of the conservation laws. More precisely, we consider the integral quantities 
formally defined by 

Ei{*) = J (M*) + hm + \), (3.12) 

= -J (um + 3/ 3 w + Ihm + \) , (3.i3) 

^3 W = jf (AW + 5/5 (*) + y/ 3 (*) + + ^) , (3.14) 



and 



f / 35 35 35 7 \ 

£ 4 W = - J r (/ 9 (*) + 7/k*) + -m*) + y/3(*) + - g-/i(*) + ye) ■ ( 3 - 15 ) 

Setting rj = 1 — 2 as usual, formal integrations by parts lead to the expressions 

E 1 (V) =E(V) = l -(\d^\ 2 + -f V 2 , (3.16) 

^(*) 4 / \9l*\ 2 -\l^\ 2 + \f (9 xV ) 2 -\fr,\ (3.17) 

4 / |d x 3 *| 2 + \ j\d 2 A 2 + \l l^| 4 + I / ^vff - I / r^ 2 *) 2 (3.18) 
z jr 4 Ji 4 jr 1 jr z Jr 

~\ [ r ] (d x r ] ) 2 + ^ [ ri 2 \d x *\ 2 + ^ [ r? 4 , 

4 JR 4 JR 1,3 JR 



and 



£4(*)4 / iW + j / l^l 2 -^ / ^) 2 -| / ^*| 2 + f / r? 2 ^) 2 

z JR 4 JR 4 JR z JR 8 JR 

?M*| 2 -f / (d^) 2 \d^\ 2 -l [ \d x *\ 2 \d 2 *\ 2 - 7 [ d 2 x n(d^,d^) 

l 4 JR z JR JR 

-7 / \ dx *\ 2 (d x *,d x *)-^ [ ndtvid^i 2 - ^ [ r^vif-^ / ^ x *| 4 

JR ^ JR 4 JR 4 JR 



35 



16 ' 11 ■ 



(3.19) 



These expressions involve only integrable integrands, and therefore provide a rigorous definition 
of the corresponding integrals. We will refer to E^l^) as the /c th -order energy. 
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When n = 2m, with m > 2, the same strategy can be applied to define the fc th -order 
momentum. We first introduce the formal linear combinations of even conservation laws 

P 2 (*) = i J (/ 4 (tf ) + |/ 2 (tt)) , (3.20) 

P 3 (v[r) = _j / (/ 6 (tf) + 5 / 4 (*) + 5/ 2 (* )) , (3.21) 
Jr j 



and 



P 4 (*) = i J (/ 8 (*) + 7/ 6 (*) + f / 4 (*) + ^/ 2 (*)) • (3.22) 

After some integrations by parts, these expressions are transformed into the well-defined quan- 
tities 

*M*) = ^ / <*^*A*> - t / (3.23) 

W = J / <i^*,^*)-| / r^d 2 *,^*) + 7 / (r^ry)^*,*), (3.24) 
z Jr z Jr 4 Jr 



and 



z Jr z Jr z Jr 4 jr 

/ v 2 (id^,d^)-^ [ (r t 3 + r t 2 + r,){id x %^). 
Jr 10 jr 



35 



(3.25) 



The case n = 2 has to be discussed separately. The invariant J 2 (\l/) is formally equal, up to 
some integration by parts, to 

This quantity is purely imaginary. Its imaginary part is equal to the momentum, i.e. 

Im(/ 2 (*)) = Px(*) = P(tf) = \ ! (id**, *). (3.26) 

2 JR 

However, the definition of the momentum raises some difficulty. As a matter of fact, the quantity 
P(^) is not well-defined for any arbitrary map ^ in the energy space A 1 (]R). We refer to [3] for 
a proof of this claim, and a discussion about the different ways to provide a rigorous definition 
of the momentum in the energy space. Notice that in our analysis of the transonic limit, we 
handle with maps ^ with small energy. In particular, we may assume that they satisfy 

E(V) < (3.27) 

so that we may lift * as 

* = gexpiip. (3.28) 
Then, we may define a so-called renormalized momentum by 

Pim= P m = l f r,d xl p (3.29) 
z Jr 
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(see [2, 3] for more details), which is also, at least formally, an invariant for the Gross-Pitaevskii 
equation, since it verifies 

pi(¥) = -i f / 2 (tt), (3.30) 
Jm. 

when ^ is sufficiently smooth and integrable at infinity. 

We will also consider the renormalized momenta p^, which are linear combinations of Pj, and 
p\. They are defined by 

P2(tf) = P 2 (V) - ^pi(*) =\ I (id^,d x ^) - - f v(id x V,y) - - A I r)d x <p, (3.31) 

z 1 JR 4 JR 4 JR 

P 3 W = P 3W + | p lW= l /" ( ^,^*)-| / ( V -l){id^,d x V) (3.32) 
z z jr z Jm. 

+ j f (V 2 + *) + t / ^x¥>, 

and 

p 4( *)^P 4W -^p lW =1 /" <^*,3 X 3 *)-| / ^ x 3 *,^) + ^ /" ^(iflg*,^*) 

o ^ JR ^ JR 4 JR 

35 /" 35 /" 

— 77 (V 3 + V 2 + rj){id x V, *) - — / fyd^, 

16 y K 16 y M 

(3.33) 

provided that the function \^ satisfies condition (3.27). As a matter of fact, the renormalized 
momenta p/,, more than the momenta P^, will be involved in the analysis of the transonic limit. 

We may summarize some of our previous discussion in 

Lemma 3.1. The junctionals Ej~, for 1 < k < 4, and P^, for 2 < k < 4, are well-defined 
and continuous on X k (R). The functionals Pk{^) are well-defined for any function ^ G A fc (M) 
which satisfies (3.27). 

Proof The proof follows from the definition of the space X 1 (R) for the functional E 1 = E. For 
the momentum p\ = p, it is proved in [2] that any function G X X (1R) such that (3.27) holds, 
verifies 

Pmin = inf > 0, 

so that, denoting *!/ = gexpitp as above, 

\vd x f\ < \rj\\gd x (p\ < — !— \rj\ |<9 X *|. 

Pmin Pmin 

Hence, the quantity rjd x (p belongs to L X (1R), so that p(^) is well-defined as well. Finally, for the 
higher order invariants, notice that, by the Sobolev embedding theorem, any function ^ G A fe (IR) 
belongs to Cq _1 (1R), so that, in particular, r\ is in H k (M). Continuity raises no difficulty. □ 

3.2 Conservation of the invariants in the spaces X fe (R) 

The purpose of this section is to provide a rigorous mathematical proof to the fact that the 
invariants are conserved along the Gross-Pitaevskii flow. As mentioned in the introduction, 
conservation of the energy E\ = E was already addressed in [20] (see also [10]). 
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Theorem 3.1 ([20, 10]). Let ^ G X 1 {R). Then, the unique solution to (GP) in 

C°(R, X 1 ^)) wii/j initial data given by Theorem 1 satisfies 

for any tel. 

Concerning the momentum, Gallo [7] established the conservation of the renormalized mo- 
mentum pi (see also [3]). 

Theorem 3.2 ([7, 3]). Let be a function in X 1 (R) which satisfies (3.27). If^(-,t) stands 
for the unique solution to (GP) in C°(R, X 1 (M)) with initial data given by Theorem 1 , then 

p{*(;t)) =p(*o), 

for any tel. 

Here, we extend the analysis to the integral quantities -Pfe(^) an d Ek{^>). 

Theorem 3.3. Let 2 < k < 4 and $ G A fc (M). T/ien, toe unique solution *(-,i) in i/ie space 
C°(M,X fc (lR)) to (GP) with initial data ^ ffwen 6y Theorem 1 satisfies 

P fe (*M)) = P fc (*o), and £*(*(■,*)) = £? fc (* ), (3-34) 

/or any i 6 R. 

Remark 3.1. Theorem 3.3 focuses on the conservation of integral quantities which play a role 
in the analysis of the transonic limit. As mentioned in the introduction, the mass m{^>) defined 
by (4) is also formally conserved. However, the quantity m(^f) is not well-defined in the energy 
space X 1 (M). A proof of its conservation along the Gross-Pitaevskii flow would first require to 
provide a precise mathematical meaning to this quantity in X 1 (M). 

Similarly, Theorem 3.3 does not address the question of the existence and conservation of higher 
order energies and momenta. A more general treatment of the inductive form of the conservation 
laws f n would be required to define properly higher order energies and momenta. However, we 
believe that such integral quantities could be well-defined in the spaces X k (R) taking linear 
combinations and integrating by parts as above, so that their conservation along the Gross- 
Pitaevskii flow would also follow from Lemma 3.2 below. 

At this stage, notice that, in view of Theorems 3.2 and 3.3, and definitions (3.31), (3.32) and 
(3.33), the quantities pk are also conserved along the Gross-Pitaevskii flow. 

Corollary 3.1. Let 2 < k < 4, and let be a function in X k (R) such that assumption (3.27) 
holds. Then, we have 

p fc (tf (■,*)) =p fc (*o), 

for any t £ 1, where ^ denotes the unique solution to (GP) in C°(R, X k (R)) with initial data 

In the proof of Theorem 3.3, we will make use of the fact that the functionals f n are conser- 
vation laws for (GP). More precisely, we have 

Lemma 3.2. Let —oo < a < b < +oo and n > 1. Consider a solution to (GP) such that 

* € C°((a,b),C n+1 (R)) nC 1 ((a,6),C"~ 1 (lR)). (3.35) 

Then, the map t ^ /„(¥(•,*)) is in C°((a, b), C 1 (R)) n C x ((a, b), C°(R)), while the function t >-> 
f n+ i(^(-,t)) belongs to C° '((a, b),C 1 (R)). Moreover, they satisfy 

dt(f n m) = id x (f n+ im - cWnW) on R x (a, 6). (3.36) 
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We will first consider the maps JF n (^) denned by (3.10), and prove 

Lemma 3.3. Let -co < a < b < +00 and n > 1. Consider a solution ^ to (GP) which 
satisfies (3.35). Then, the map t ^ T n (^(-,t)) is in C°((a,b),C 1 (R)) nC 1 ({a,b),C°(R)), while 
the function 1 1— > !F n+ i{^> '(■ ,t)) belongs to C° '((a, 6),C 1 (M)). Moreover, they satisfy 

n-l 

dt(F n (*)) = £F„(tf) - i|*| 2 ^n(*) + tftc* ^^(*)^ n _j(*) + id^Fn+im). (3.37) 

j'=i 

Lemma 3.2 is then a direct consequence of Lemma 3.3. 

Proof of Lemma 3.2. Notice first that, in view of assumption (3.35) and formulae (3.2) and 
(3.10), the maps and belong to C°((a, b), C 1 (R)) for any 1 < j < n + 1, and the 

functionals f n (^) an d F n (^) are also in C x ((a, 6), C°(R)). Therefore, in view of (3.8), we can 
write 

$(/„(*)) = dt*r n (*) + *dt(r n {*)), 

so that, by (GP) and (3.37), 

n— 1 

(/"(*)) = *(- 3^n(* ) + Wj^f j ($)f n _ J ($) +^9 X (^„ +1 (*))). 

J'=l 

In view of (3.10), we are led to 

^ (/„(*)) = i( " ^*^„(*) + d x ^Tn + im - dJ£d x {F n (y)) +Vd x (F n+1 (*))), 

which completes the proof of (3.36), invoking definition (3.8). □ 

We now provide the proof of Lemma 3.3. 
Proof of Lemma 3.3. The proof is by induction on n € N*. For n = 1, it follows from (GP) that 

flt(.Fi(¥)) = -iflt* = 1(- d 2 * - * + |^| 2 ^) = tfi(tf) - i|tf| 2 .Fi(¥) +^ X (JP 2 (*)), (3.38) 

so that (3.37) holds for n = 1. We now turn to the case n = N + 1, assuming that the conclusion 
of Lemma 3.2 holds for any 1 < n < N . Notice first that, in view of assumption (3.35) and 
formulae (3.2) and (3.10), the maps are in C°((a, b), C X (R)) for any 1 < j < N + 2, while 

the functional .F/v+i^) also belongs to C 1 ((a, b),C°(M)). Therefore, in view of (3.10), we can 
write 

JV-l 
3=1 

N-l 
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Invoking the inductive assumption combined with (GP), we are led to 

d t {f N+1 (v)) =i(ji - \*\ 2 )d x (F N (t>)) - (vajz + *d x *)r N (*) 

N-l 

+ *(i - 1 if) e rm?N-jW 

3=1 

N-l 

3=1 
N-l 

3=1 

N-l j-1 N-j-1 . 

+$W^(^. J ($)^^($)f rfc ($) + f ; (f) £ 

7 = 1 fc = l fc = l ' 



(3.39) 



In view of (3.9), we first have 

JV-1 



J'=l 

AT 

=*E +^(«)9x(^ +H w)), 



(3.40) 



whereas, by formula (3.10), 

N-l 

(i - |*| 2 ) (d*{r N (*)) +*Y. rjWFN-jm) = (i - |*| 2 ).w W, (3.4i) 

3=1 

and 

JV-1 j'-l N-j-1 

-VdxVFNi*) + Vd x V J2 (-^v-i W E ^W-^-fc W + ^ (*) E ^fc(*)^v-j-fc(*)) 

i=l fc=l k=l 

N N-l 

(3.42) 

Hence, we deduce from (3.39), (3.40), (3.41) and (3.42) that 

N 

E 

(3-43) 



dt(r N +i(*)) =i[ (i-|*l 2 )-WW + d^Y^^N+i-jm 

3=1 
N 

+d x (d x {T N+1 E FjWfN+i-j (*)) 



j'=i 



In view of (3.10), the second line in (3.43) is equal to d K (j r N+2 ( i ^)) , so that (3.37) holds for 
n = N + 1. This completes the proof of Lemma 3.3 by induction. □ 
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We finally turn to the proof of Theorem 3.3 

Proof of Theorem 3. 3. We first assume that in addition £ X 9 (R) . In this situation, the maps 
t i ^ E k ($>(-,t)) and t ^ P k (^(-,t)) are in C^IR,^), while by the Sobolev embedding theorem, 
the map t i-> is in C 1 (M,C 8 (M)) and C°(R, C 10 (R)). Hence, in view of Lemma 3.2, 

9t (/„(*)) = iftc (/„ + i(*) - S^C*)) on R, (3.44) 

for any 1 < n < 9. 

Now consider, for instance, the map t ^ £7 2 (\l/(-, i)). In view of (3.13), its derivative is, at 
least formally, given by 

^(Eh^t))) =-J (dtuw+mhm + 

so that by (3.44), we formally have 

A(^(*(.,t))) =-i|^^(f) + 3/ 4 (*) + ^(*)-W(f5W + 3f2W + ^iW)) =0, 

i.e. the quantity E^C 1 ^) is formally conserved by (GP). In particular, the proof of the conser- 
vation of £2 along the Gross-Pitaevskii flow reduces to drop some integrability difficulties in 
the above formal argument. Therefore, given any R > 1, we introduce some cut-off function 
X G C°°(R, [0, 1]) such that 

X = 1 on (-1, 1), and x = on IR \ (-2, 2), (3.45) 

and denote 

Xfl(x) = xQ,V X GR. (3.46) 
Since the map 1 1-> belongs to C X (1R, X 9 (M)), we then have 

A(E 2 (f(-,t))) = / dt(e 2 (*(;t))) = lim / X fl(x)^(e 2 (*(x,t)))dx, (3.47) 
where we let 

E2W = [ e 2 (Y>). 
Jr 

We now make use of formal relation (3.13) to compute 
/ Xfl(x)$(e 2 (*(x,t)))<k = - / X fl(^/5(*) + W 3 (*) + |^/i(*))+ / d xXR Qi(*,$*), 

JR V ^ 7 JR 

where, using definitions (3.1), (3.5), (3.7) and (3.17), and the Sobolev embedding theorem, the 
function Qi(^,dt^) tends to at ±00. Invoking (3.44) and integrating by parts once more, we 
are led to 

/ Xfl(x)9t(e 2 (tf(M)))<k= / d xX RQ2(V,dt*), (3.48) 

JR JR 

where 

Q 2 (* , = ) + */e(*) + 3i/ 4 (*) + |*/ 2 (*) - id x *(V 5 (*) + &F 2 («) + , 

also tends to at ±00. Finally, notice that when 

/(x) — > 0, as |x| — > +00, 
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we have 



J^di XR (x)f(x)dx= ^rr(^ 2 a£x(x)/(ife)dx + j \{ X {x)f{Rx)d^j -0, asi?^+oo, 
so that, in view of (3.47) and (3.48), we obtain at the limit R — > +oo, 

5(a(»(-,*)))-o, 

which gives (3.34) for the quantity 

Using formal identities (3.14), (3.15), (3.20), (3.21) and (3.22), the proofs are identical for 
the functional E3, E4, P2, P3 and P4, so that we omit them. 

In the general case where we only have ^0 £ X fc (R), we first approximate ^0 by a sequence 
of functions ip n in X 9 (1R) for the X fc -distance (see e.g. [9]), and then use the continuity of the 
flow map $0 l— * ^(',T) in X fe (M) for any fixed T, and the continuity of the functionals E k and 
Pk with respect to the X fc -distance. □ 



4 Invariants in the transonic limit 

In this section, we analyse the expressions of the invariant quantities introduced in the previous 
section in the slow variables. Therefore, we introduce the quantities £k(N £ , e ) and Vk(N £ , Q £ ) 
defined by 

£ 2k+l 



E k (V) = —£ k (N £ ,@ e ), (4.1) 



and 



£ 2k+l 



We also set 



PkW =^V k (N £ ,G £ ). (4.2) 

m £ = 1 - — N £ . 
6 

We now derive the precise expansions of £ k and Vk and stress the relationship with the corre- 
sponding (KdV) invariants. 

4.1 Formulae of the invariants in the rescaled variables 

For the fc t/l -energies defined by (3.17), (3.18) and (3.19), a direct computation provides, in view 
of definitions (14) of N £ and Q £ , 

Lemma 4.1. Let 2 < k < 4 and e > 0. Given any function ^ in X k (M.) which satisfies (3.27), 
and denoting N £ and £ , the functions defined by (14), we have 

£ 2 (jv„ e.) = i J ((om) 2 + m , (die, - ^a.jv.a.e,) 2 - ijvf - ^-n, (d l: e s f\ 

with 

R2W ,e e)s W£>! ,4.4) 
m F 
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1 f ( ^3n *\«r^3 e 2 d x N £ d 2 x Q £ e 2 d 2 N £ d x Q £ e\d x N £ ) 2 d x Q £ ^ 



4m £ 4m e 48m| 



+ (^iv £ ) 2 - -(2d x N £ d x e £ d 2 e £ + N £ {d 2 x e £ ) 2 + N £ {d x N £ ) 2 ) + _((^e £ ) 4 

+ m 2 {d x Q £ f + A^ 4 )) + ^ jf -L " y 4 d*Md*®c) 2 + ^d x Q £ d 2 x Q £ 

(4.5) 



where 



K 3 (N £ ,e £ ) = lN 2 {d 2 x Q £ f - b -{d x N £ f{d x Q £ f - 5N £ d 2 x N £ {d x e £ f - ^N*{d x G £ ) 2 
~ r 8 N ^e) 4 - ^(0 X@£ ) 2 - ^-N £ {8lN £ ) 2 + ^-N 2 {d x N £ f + t {^N 2 {d x Q £ f 
~ ^Md x N £ ) 2 (0 x e £ ) 2 - ^(0 X N £ ) 4 - ^N £ (d x N £ ) 2 0lN £ - ^N £ {d x N £ f). 

(4.6) 

and 

- £ ^(d x e £ ) 2 d 2 x e £ - ^(d,N £ ) 2 d 2 x e £ - -^d x N £ d 2 N £ d x e £ + ^J^M^f 

~ T^S dxN ^ dxQe ) 2 + < W -1{ N °& N <) 2 + m £ N £ (dl@ £ ) 2 

+ 2m £ d 2 x N £ d x <d £ dl<d £ ) + ^N 2 (d x N £ ) 2 + m £ (d x e £ ) 2 (d 2 e £ ) 2 + m £ d x N £ d x & £ d 2 x e £ 
+ m £ N 2 {d 2 x e £ ) 2 + {d x N £ ) 2 {d x e £ ) 2 ) - ^(at £ 5 + siv^e,) 4 + ioiv|(a x e e ) 2 ) 

+ 15 + + f + 3^^ + 4^(« 2 

(4.7) 
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where 



35 , t q,„ ^ N 4 , 7 /r , „ A 6 35 ,w^^\2/ ^\2 



n A {N £ ,Q £ ) = —m £ N 2 £ {d x Q £ y + — m £ (d x Q £ y - -m £ N £ (d x Q £ y(d x e £ ) 

- ^N%N £ d x @ £ dl@ £ + f 2 dlN £ (d x e £ f - l ^N £ (d x N £ )\d x e £ f - f A {d x N £ )\dl@ £ f 
+ ^(d 2 x N £ ) 2 (d x & £ ) 2 + ^-d x N £ d 2 x N £ d x e £ d 2 x e £ + ^N £ d 2 x N £ (d 2 x e £ ) 2 + ^N £ d 2 x N £ d x @ £ d 3 x @ £ 

7, " £ iv £ (^e £ ) 6 - §A^iv £ (a,e £ ) 4 - 5j(a,iv £ ) 2 (a,e £ ) 4 + ^-iv £ 2 (a,iv £ ) 2 (a,e £ ) 2 



46656 tv c/ 54 a 2 c v " ty 432 v * c/ v * 72n\ 
+ 3^(^) 2 (^) 2 " -^Md x N £ ) 2 (d 2 e £ ) 2 - ^{d x N £ fd 2 x N £ {d x Q £ f 

^-N £ (d 2 N £ ) 2 (d x e £ ) 2 - ^-N £ d x N £ d 2 x N £ O x @ £ d 2 x @ £ - ^N £ {d x N £ f) 



M£l^ N ^^ 2 - g^^(^) 2 (^e e ) 4 + §^(^) 2 ^(^e £ ) 2 

& 2 (^) 4 + ^ w.) a (w a ) + 4(^(^) 4 ^) 2 



497 ,„ „ n6 \ 7e 8 



(4.8) 



Similarly, for the fc^-renormalized momenta, we compute 

Lemma 4.2. Let 2 < k < 4 one? £ > 0. Given any function ^ in X fc (R) which satisfies (3.27), 
and denoting N £ and £ , the functions defined by (14), we have 

P 2 (iv £ ,e £ )^^^(^ £ ^e £ -^(a,e £ ) 3 -^ £ 2 ^e £ ) r2 (jv e> e e ), (4.9) 

r 2 (iV £ ,e £ )^^^^, (4.10) 
m £ 

p 3 (iv e ,e £ ) =^= jf (^iv £ ^e £ - ^(^e^e,) 2 + 2N £ d x N £ d 2 x e £ + (sjv,) 2 ^) 
+ 4(m £ iv £( ^e £) 3 + a^ 0£ )) + / /3 (iv £ ,e £ ), 

(4.11) 

with 

r 3 (iv £ ,G £ ) = ^iY £ ^e £ (^e £ ) 2 + h x N £ {d x ® £ ) 2 d 2 x ® £ - ^{d x e £ ) 5 + ^-7v £ (^7v £ ) 2 ^e £ 



~ ^N e d 2 x N £ d 2 x e £ - ^- £ {dlN £ fd x Q £ - ^{d x N £ f(d x ® £ f - -^{d x N £ fd 2 x Q £ 



(4.12) 
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and 

v 4 (n £ , 9 £ ) = -LJ^ (dlN £ d A x ® £ - ^ (2N e d 2 x N £ dl® e + (d 2 x N £ ) 2 d x o £ + m £ d x Q £ {dle £ f 

(N 2 £ d x N £ d 2 x @ £ + N £ (d x N £ ) 2 d x e £ + d x N £ (d x e £ ) 2 d 2 e £ + m £ N £ d x e £ {d 2 x o £ ) 2 ) 



35 
+ T2 



- Y^((^e £ ) 5 + iom £ N 2 {d x e £ ) 3 + 5jv^a x e e )) + ^=^r 4 (iv £ ,e £ ), 

(4.13) 

with 



r 4 (iv £ ,e £ ) = 7(— a^at,) 2 ^ + — A^iv £ d 2 iv £ <9 2 e £ - ^-N*(d x N e ) 2 d x e e 
+ -^iVAe^e,) 2 - ^A^iv £ (^e £ ) 2 d 2 e £ + ^ 2 iv £ (^e £ ) 2 ^e £ - ^(d x N £ ) 2 (d x & £ ) 

+ 2l6^ae £ ) 5 - y 2 d x N £ {d 2 x Q £ f - ^{d x Q £ f{d 2 x Q £ f - ^- e {d x N £ fd 2 x Q £ 

+ 7L~S dxN ^^? - ^d x N £ d x N £ d x e £ - -±-(d*N £ ) 2 d x & £ + ^(d x N £ fd 2 x ® £ ) 

+ -(^(dM 4 d x e £ - -^N 2 (d x e £ ) 5 - ^(d x e £ ) 7 - -^(d 2 N £ ) 2 (d x e £ ) 3 

4V432m2 v £j x 6 1296 eV ; 1296 v ; 36m £ v x £ > K £ > 

35 a at f>2*r fa o ^a2^ 35 /a Ar \2« ^ /o ^ \2 , 7 Ar /a2 Ar \2 Q 2 
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■a x jv e ^iv e (a x e e ) 2 ^e e - — (5xJv e ) ^e e (a x e e ) 2 + ^d x N £ (d 2 x N £ ) 2 d 2 x @ £ 
7 (d x N £ ) 2 d 2 x N £ dle £ + ^(a 2 iv £ ) 3 ^e £ - ^(9 x n £ )%@ £ - -^d 2 N £ (d x e £ ) 5 ) 



ml 

+ i {^ x N £ fd 2 x N £ d 2 x e £ - -t-{d x N £ ) 2 {d x ® £ f -^(d x N £ f( dx e £ ) 2 d 2 x e 



48\ml y x £J 7 2m £ v x £ > v x £ > Qm? £ 

~ iL £ (°° N <) 2 % N <fre<) 3 + ^ (^ £ ) 2 (a 2 iv £ ) 2 ^e £ ) + ^(^(^) 4 (^e £ ) 3 
+ ^(^ £ )^ 2 e £ -i^(^ £ )^e £ ). 

(4.14) 



4.2 Relating the (GP) invariants to the (KdV) invariants 

Recall that the Korteweg-de Vries equation is integrable, and admits an infinite number of 
invariants (see [8]). The first four invariants for (KdV) are given by 

< d » 4 / v 2 , (4.15) 

£f d »4 / &V) 2 - 1 - [ v\ (4.16) 

E? dV (v) ^ / (a 2 ,) 2 - \ [ v(d x v) 2 + ^fv\ (4.17) 

J 1R •/ 1R t/IK 



and 



Notice that the invariants are bounded in terms of the ff fc -norm, since we have 

\E% dv {v)\ < K(\\v\\ Hk - 1{R) )\\v\\ 2 Hk{R) , (4.19) 
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where ^(IMI/f^-^R)) is some constant depending only on the H k 1 -norm of v. 

Another important observation concerning the (KdV) invariants is that, given any function 
v G H k (W), the if fc -norm of v is controlled by the first fc^-invariants of (KdV). This claim is 
straightforward for k = 0, whereas for k > 1, we have 

Lemma 4.3. Let 1 < k < 3 be given. Given any function v £ H k (W), there exists some positive 
constant K = K (\\v\\ H k-i ^) , depending only on the H k ~ l -norm of v, such that 

\\d k x vf L2(m < K(\\v\\ 2 Hk . 1{m + \E« dV (v)\). (4.20) 

Proof. For k = 2 and k = 3, the proof of (4.20) is a direct application of the Sobolev embedding 
theorem to formulae (4.17) and (4.18). For k = 1, we have in view of (4.16), 

\\d x v\\ 2 L 2 m < l\\v\\h m \\d x v\\ L 2 {R) +2\Ef dV (v)\, 

so that by the inequality 2ab < a 2 + b 2 , 

\\d x v\\h m < l\Ml Hm +4|^ d »l < K {\Hh m + 

where K = max {|||i;||| 2 ^,4}. □ 

We complete the subsection showing that the (KdV) invariant Ejf^Y is related to the (GP) 
invariant quantities ± \[2Vk- For that purpose, assume that 

N £ -► V in if^K), and c^Be -► ^e in L 2 (M), as e -► 0. 

For k = 1, we notice in view of expansions (28) and (32), that 

gi(JV e , 6 e ) ± V2Pi(V £ , 8 e ) -» ±J (N ± d x O ) 2 = E« dV ( ^ ± 2 dx@0 ) , ase^O. (4.21) 

Similarly, it follows from Lemmas 4.1 and 4.2 that 

Proposition 4.1. Let 1 < k < 4 and ^ m X k (M.) which satisfies (3.27). Denoting N £ and Q £ 
the variables defined by (14), and assuming that 

N £ 7V in H k (R), and c\0 £ -»• 9 X O in ^ fe_1 (M), as e ->■ 0, 

we /iaue 

£ fc (iv £ ,e £ ) ± V2V k (N £ ,@ £ ) - Etf( No± ^ xeo ), as e - 0. 

Remark 4.1. We believe that Proposition 4.1 might be extended to higher order (GP) and 
(KdV) invariants, provided one was first able to compute some expressions for them. 

Proof. Combining the expansions of Lemmas 4.1 and 4.2 with (4.16), (4.17) and (4.18), and 
using the Sobolev embedding theorem, the proof reduces to a direct computation similar to the 
proof of (4.21). □ 
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4.3 //^-estimates for N £ and d x Q £ 

In the same spirit as Lemma 4.3 which allows to bound the i/ fc -norms by the (KdV) invariants, 
we next show that the i? fc -norms of N £ and d x Q £ are controlled by the quantities Sk{N £ , Q £ ) in 
the limit e —>■ 0. More precisely, we have 

Lemma 4.4. Let 1 < k < 4 be given, and assume that there exists some positive constant A 
such that 

Sj{N e ,Q e )<A, (4.22) 

for any 1 < j < k. Then, there exists some positive numbers ea and Ka, possibly depending on 
A, such that 

ll^ll^-i(R) + e||^JV e || £2(R) + \\d x Q e \\ Hk -, {K) < K A , (4.23) 

for any < e < e A . 

Remark 4.2. We again believe that Lemma 4.4 might be extended to higher order (GP) and 
(KdV) invariants, which will provide bounds for higher Sobolev norms of N £ and d x Q £ . 

Proof. We split the proof in four steps according to the value of k. 
Step 1. k = 1. 

In view of (28), assumption (4.22) may be written as 

>2 , „2 , £ 2 «W 2 



/ 



(m £ (d x Q £ f + iV £ 2 + £ ^ ) < 8A, (4.24) 



so that (4.23) follows once some lower and upper uniform bounds on m £ are established. Indeed, 
if we can choose ea so that 

2 < me < 2, (4.25) 

for any < e < e A , then (4.23) follows from (4.24) with Ka = 24\/A Hence, the proof reduces 
to show (4.25) for some suitable choice of ea- 

In order to prove (4.25), we apply the Holder inequality and assumption (4.22) to obtain 

\V(x)-V(xo)\ < ^x-xq^E^ < ^(iV £ ,e £ )3 < ^_ e §, 
for any point xq £ R and any xq — 1 < x < xq + 1, so that 



|l-|*(x )||-^et<|l-|*(x)||. 

Setting ea = (16A) _ 3 5 and assuming by contradiction that (4.25) does not hold at the point xo, 
we obtain that 

|1 - |*(x )|| = |1 - y/m e (x )\ > 1 - > ^ = ^e\ > ^e 3 , (4.26) 
for any < e < e A , so that 

(|l-|*(xo)||-^S) <j x( i {l-^ix^fdxK^N^)^^, 
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and 

|l - |*(zo)|| < >/Ae* < _ 

It follows that 

9 25 
< m £ (x ) = \^{x )\ 2 < —, 

which gives a contradiction with the fact that (4.25) does not hold at the point xq. This 
completes the proof of (4.25), and of Step 1. 

Step 2. k = 2. 

Notice first that in view of the inductive nature of assumption (4.22), and of Step 1, we 
have already established (4.23) for k = 1. Combining this estimate with the Sobolev embedding 
theorem, bounds (4.25) and formulae (4.3) and (4.4), assumption (4.22) may be written as 

< K(l + \\N £ \\ Him (\\N £ \\l 2m + Re e ||! 2(R) + £ 2 ||^iV £ ||| 2(R) )) < K A (l + \\N £ \\ m{m 

(4.27) 

This first gives that 

/ {d x N £ f<K A , (4.28) 
Jr 

so that by (4.27) and the Sobolev embedding theorem, 

\\dl@ £ \\ L 2 m <K(i + e 2 \\d x N £ d x e £ \\ L 2 m ) <K A (i + £ 2 ||^e £ ||^i (ffi) ). 

Hence, we obtain 

(d 2 x ® £ ) 2 < K A , 



J 

JR 



setting ea sufficiently small. In view of (4.27), (4.28) and the Sobolev embedding theorem, it 
follows that 

£\\dl N e\\tf(R) < ^(l + e||^0 £ |li4 (R) +e 3 ||^A r £ ||//i(R)||^iV £ || i 2 {]R) ) < K A (l + e 3 \\d x N e \\ H i {R) ^, 
which completes the proof of (4.23) choosing ea sufficiently small. 
Step 3. k = 3. 

Notice again that in view of the inductive nature of assumption (4.22), and of Step 2, we 
have already established (4.23) for k = 2. Combining this estimate with the Sobolev embedding 
theorem, bounds (4.25) and formulae (4.5) and (4.6), assumption (4.22) may be written as 

f ((#n\ 2 +(&* g2 (a ft ^ e 2 d x N £ d 2 x Q £ E 2 d 2 N £ d x @ £ E\d x N £ fd x Q £ ^ 
LV X ' + \ x 6 72 [dxUe) 4m £ Am £ 48m? ) 

+ e 2 (dlN £ - £ l dx N £ {d x e £ f + ^d x e £ d 2 e £ + ^- dx N £ d 2 x N £ + J^( dx N £ f) 2 

<K A (l + \\d x N £ \\ mm + \\d 2 x e £ \\ mm +e\\a%N e \\ m 



Invoking once again estimates (4.23) (for k = 2) and (4.25) to bound the remainder terms in 
the above integral, we are led to 



{d 2 x N £ f + (dl® £ ) 2 + E 2 {dlN £ ) 2 ) <k a (i + \\d x N £ \\ mm + \\d 2 x @ £ \\ mm + E\\dlN £ \\ L 2 { ^, 

which provides the proof of Step 3. 
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Step 4. k = 4. 

Notice once last time that, in view of the inductive nature of assumption (4.22), and of Step 
2, we have already established (4.23) for k = 3. Combining this estimate with the Sobolev 
embedding theorem, bounds (4.25) and formulae (4.7) and (4.8), assumption (4.22) may be 
written as 



(d 3 x N £ ) 2 + (d x O £ - ^d 2 N £ d 2 x Q £ - ^d x N £ d x @ £ - ^-dlN £ d x Q £ 



~ ^{d^efdlQe - ^{d,N £ fdlQ £ - ^d x N £ d 2 x N £ d x ® £ + i^d x N £ {d x Q £ f 
~ 7^3 {d X N £ fd x Q £ ) 2 + e 2 jf (d x N £ + ^(d 2 x O £ ) 2 + 2 -m £ d x ® £ dl® £ + J^-d x N £ d x N £ 

+ -& £ ^ Ne ^ ~ i^7v £ (3 x £ ) 2 - £ ^d x N £ d x e £ d 2 x e £ - S^m £ (d x e £ ) 4 



£ 



(d x N £ ) 2 (d x Q £ ) 2 + ^ (d x N £ ) 2 dlN £ + ^ (9,iV £ ) 4 ) 2 



144m o 

<K A (l + ||^7V £ ||^(ffi) + ll^© £ |ki 
so that we similarly obtain 

/ M (a 3 iv £ ) 2 +a 4 e £ ) 2 )<^, 

then, combining with the Sobolev embedding theorem, we also have 

4d x N £ \\ L 2 {R) < K A . 

This completes the proofs of Step 4 and Lemma 4.4. □ 

An important consequence of Lemma 4.4 which refines the result of Lemma 4.1 is 

Proposition 4.2. Let 1 < k < 3. Given some positive constant A, consider some functions N £ 
and d x @ £ which satisfy (4.22) for any 1 < j < k + 1. Then, there exists some positive numbers 
ea and Ka, possibly depending on A, such that 

£ k {N £ , 9 e ) ± V2V k (Ns, 6 e ) - Egy ( - ± / xee ) I < ^ 2 - (4.29) 
/or any < e < . 

Remark 4.3. Similarly to Lemma 4.4, we believe that Proposition 4.2 might be extended to 
higher order (GP) and (KdV) invariants. 

Proof. Let k = 1. In view of (28), (32) and (4.15), we have 

£l (N £ ,@ £ ) ± v^w,e £ ) - ^(^MMi) = * _ i^^e.) 2 ). 

Inequality (4.29) follows for sufficiently small, invoking (4.23) (for k = 2) and (4.25). 
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For k = 2, we deduce from (4.3), (4.4) and (4.16) that 

£ 2 (n £ , e e ) ± V2v 2 (n £ , e e ) - < dy ( iYe± 2 ^ 9£ ) 
= j jf (^(a*e e ) 2 - \d x N £ d x ® £ dl® £ - l -N £ {dl® £ f ± ^iv £ (^e £ ) 3 - 



4m £ 



T 4m £ + 36m £ 2m £ \ x 6 6 + 12m £ [ ° x e) ) )' 

so that (4.29) follows again from (4.23) (for k = 3) and (4.25). 

Similarly, the proof of (4.29) for k = 3 reduces to estimate the remainder terms in (4.5) and 
(4.6) using (4.23) (for k = 4) and (4.25). □ 

5 Time-independent estimates 

In this section, we use the above conservation laws to derive time-independent estimates of the 
functions U £ and V £ , together with the consistency of the solutions to (GP) with the (KdV) 
equation in the limit e — ► 0. This yields the proofs of Proposition 2 and Theorem 5. 

5.1 Proof of Proposition 2 

Given any functions N £ and @ £ such that (15) holds, it follows from the formulae of Lemmas 
4.1 and 4.2 that there exists some positive constant Aq, which does not depend on e, such that 

£ k (N^Q° £ )<A , (5.1) 

for any 1 < k < 4. In view of Theorem 3.3 and definition (4.1), we deduce that the solution 
(iV £ (-,r),e £ (-,r)) to system (23)-(24) with initial datum (iV £ ,e°) satisfies 

£ k (N £ (;T),@ £ (;T)) <A , 

for any time r£l. In particular, inequality (36) is a direct consequence of Lemma 4.4, whereas 
in view of Proposition 4.2, we have 

£ k (N £ (;r),e £ (.,r)) ± y/2P k (N e (; r), £ (, r)) - [ NeM ± ^^ ) | < K M e\ 

for any time rgl. Using again the conservation of E k and p k provided by Theorem 3.3 and 
Corollary 3.1, and definitions (4.1) and (4.2), we are led to 

£ k (Nl e°) ± V2v k (Nl e°) - Et d r (Mill^Mil)) I < KmE \ 

Invoking (5.1), we apply once more Proposition 4.2 to obtain 

KdV ( N° ± d x eh „ dV ( N £ (• , r) ± d x £ (• , r ) 



rpi^av i s e l 171/ 

I rj ' ~~ fc ~ 3 



< K Ao e 2 . (5.2) 



For k = 1, we then deduce from (4.17) that 

||iV £ (-,r)±^e £ (-,r)|| L2(ffi) < ||Ar°±^e°|| L2(M) +K Ao£ , (5.3) 
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so that in particular, we have by (15) for e sufficiently small, 

\\N e (-,r) ± d x Q £ (-,T)\\ L 2 m < K Ao , 

where Ka denotes some further constant depending only on A$. Hence, for k = 2, we may 
write using (4.20), that 

||^ £ (-^)±^e £ (.^)||^ 2 < ^ (|^f^( ^ £( " ; ' ) ^ e£( '^ ) )| + ||Ar £ (.^)±^e £ (-,x)|| z , a ), 

so that by (4.19), (5.2) and (5.3), 

\\d x N £ (;T)±d 2 x e £ (;T)\\ L 2 (R) < K A() (||iV £ ± ^ 1 1 ^ (R) + £ ) . 

Using repetitively this argument to estimate the L 2 -norms of the functions d 2 N £ (-, t)±c^O £ (-, t) 
and d x N £ (-,r) ± d x Q £ (-,r), we are led to (37), which completes the proof of Proposition 2. 

5.2 Proof of Theorem 5 

Theorem 5 is a consequence of Proposition 2. Applying estimates (37) to the right-hand side of 
(25), together with the Sobolev embedding theorem, we obtain estimate (22). 

6 Energy methods 

This section is devoted to the proofs of Theorems 3 and 4, which both rely on applying standard 
energy methods to equations (25) and (26). 

6.1 Proof of Theorem 4 

In order to estimate the L 2 -norm of V £ (-,t), we multiply equation (26) by V £ (-,t) and integrate 
by parts. In order to simplify the presentation, we recast equation (26) as 



d T V £ + ^8 X V £ = \d x {V 2 ) + d x f £ + e 2 R £ , (6.1) 
e z 2 



where 



f £ = d 2 x N £ - l -U 2 £ - ±U e V e , 



and R £ is defined in (27). We are led to 

d T ( f V £ (-,t) 2 ) =-2 [ f £ d x V £ (-,r) + 2e 2 [ R £ (-,t)V £ (-,t). 
We now integrate with respect to the time variable to obtain 

/ (V £ (-,t)) 2 = [ (V £ °) 2 -2 [ T [ f £ d x V £ + 2e 2 f f R £ V £ . (6.2) 

JR JR JO JR JO JR 

Combining inequalities (36) with definition (27) and bound (30) and using the Sobolev embed- 
ding theorem, we next have 

ll^(-,r)|| H 3 W + \\V £ (;T)\\ H s m + ||i? £ (-,r)|| L2(K) + ||/ e (-,T)|| H iOR) < K, (6.3) 
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for any r£l and some positive constant K not depending on e. In particular, 



2e z 



f'T P f'T 

/ / ReVe <CS 2 / \\V e (;8 

Jo Jr Jo 



Il 2 ( 



\ds 



(6.4) 



where C = C(K) does not depend on e. In order to bound the second term in the right-hand 
side of (6.2), we replace the quantity d x V £ in (6.2) according to (6.1), so that 

T / fed x v £ = ^ f T [ fj- d T V £ + \d x {V 2 ) + d x f £ + e 2 R £ ) 
Jo Jr 8 Jo Jr v z j 

= Ji + J 2 + J 3 + J±- 

We bound each of the terms separately. First note that the integrand being a differential, 
J3 = 0. Next, it follows from (6.3) that 



|J 4 | < Ce 4 T. 



Concerning J2, we have 

.2 f'T 



U2 



£ 

16 



r [ f £ d x (v £ 2 ) = ^ r f d x f £ v e 

Jo JR 10 JO 



\Ve(-,s)\\ L 2( R )ds 



For Ji, we perform an integration by parts with respect to the time variable, so that 



Note that by (23), (25), (36) and (6.3), 



/ 

Jr 



feV £ 



(6.5) 



(6.6) 



(6.7) 



J 



d T f £ = d 2 x d T N £ - hj £ d T u £ - \d T u e v £ - hj e d T v e 
4 



.^dlv £ -\u £ d T v £ + o(i) 



uniformly in L 2 (R), so that 

.2 fT 



fl 

Jo Jr 



drfeV £ 



< 



^ f j U £ d T (V £ ) 2 +Ce 2 [ T \\V £ (;s)\\ LHR) ds 
Jo Jr Jo 



£ 

48 



A further integration by parts in time leads to 

-2 rr r r-1 rr 



£ 
48 



[ T [ U £ d T (V £ ) 2 = - £ - f I (d T U £ )V 2 + £ - [ u £ v 2 
Jo Jr 48 Jo Jr 48 Ljr 



(6.8) 



(6.9) 



and since d T U £ is uniformly bounded in L 2 (IR) by (25), (36) and (6.3), we obtain, combining 
(6.7), (6.8) and (6.9), 



Ji<Ce 2 (\\V £ (-,0)\\ L 2 m + \\V £ (-,T)\\ LHR) + £ \\Ve(;s)\ 
Finally, combining (6.2), (6.4), (6.5), (6.6) and (6.10), we obtain 
\\V £ (;r)\\ 2 L2 < \\V £ (;0)\\ 2 L2 + Ce 2 (e 2 r+\\V £ (;0)\\ L 2 + \\V £ (-,t)\\ l i + 



(6.10) 



\V £ (;s)\\ L 2d S 



The proof of Theorem 4 then follows by the Gronwall lemma. 
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6.2 Proof of Theorem 3 

We first recall the equation (25) satisfied by U £ , namely 

d T u £ + dlu £ + u £ d x u £ = -dlv £ + l -d x (u £ v £ + ^f)- e 2 R £ , 

and take the difference with the (KdV) equation 

d T Ne + dlM £ + KdxAfe = 0, 

so that Z £ = U £ — N £ satisfies the equation 

d T z £ + dlz £ + z £ d x u £ + M £ d x z £ = -& x v £ + hs x (u £ v £ + ^f)- e 2 R £ . (6.11) 

We multiply (6.11) by Z £ , integrate on R and perform an integration by parts to obtain 

^r||^£|||2( K ) < K(^\\d x U £ \\ L oo(^ + 1 1 9 x J\f £ 1 1 £,oo QfJ) ) 1 1 Z £ ||^2 

+-^l|-^e||L 2 (R) (||K||h 3 (R) + l|K||L2(R)(||^e||H 1 (R) + II^eIIh^R)) + £2 \\ R e\\L 2 (R)Y 

Using bounds (6.3) for U E , V £ and R £ , and the bound of M £ in i? 3 (R) which follows from the 
integrability theory of (KdV), we are led to 

di-H^ell^n) < K||Z e ||^ 2 ( R ) + i^||^ £ ||i2( R )(e 2 + ||K||h3(R))- 

Finally, we invoke Proposition 2 to assert 

d-rll^H^^) < i^||Z e || 2 2 ( R ) + -ftT||Z e || L 2( R ) (e + ||V r £ (-,0)|| i? 3( M )), 
so that by the Gronwall lemma, 

WZeMWvw < ||Z £ (-,0)|| L 2 (R) + K( e + ||l/ £ (-,0)||^3 (M) )exp(^T). (6.12) 
On the other hand, at time r = 0, since Af £ (-, 0) = N £ (-, 0), we have 

\\Z £ (;0)\\ L2m = \\U £ (-,0) -A4(-,0)|| L 2 (M) = ||y £ (-,0)|| L2(M) , (6.13) 
whereas at positive time, by definition of V s , we have 

\\N £ (;T)-M £ (;T)\\ L 2 m < \\Z £ (-,T)\\ L 2 (R) + \\V £ (;T)\\ L 2 (R) 



< \\Z £ (;T)\\ L 2 iR) + ||K(-,0)|| L 2 (M) + Ke 2 \r\, 



(6.14) 



where we have used Theorem 4. The conclusion for N £ — J\f £ then follows from (6.12), (6.13) 
and (6.14). The proof is similar for d x Q £ — M. £ considering the function Y £ = U £ — M. £ instead 
of Z £ , so that we omit it. 
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